Stratified horizontal flow in vertically vibrated granular layers 
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A layer of granular material on a vertically vibrating sawtooth-shaped base exhibits horizontal flow 
whose speed and direction depend on the parameters specifying the system in a complex manner. 
Discrete-particle simulations reveal that the induced flow rate varies with height within the granular 
layer and oppositely directed flows can occur at different levels. The behavior of the overall flow is 
readily understood once this novel feature is taken into account. 
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Understanding the static and dynamic properties of 
granular materials is a challenging task jl| and many as- 
pects of their behavior are unexpected. An example of 
this kind of behavior is provided by a granular layer that 
is vibrated vertically by a base whose surface profile has 
a sawtooth form; in this particular case the material is 
confined to the annular region between two upright cylin- 
ders and is thus practically two-dimensional. What is 
observed j^, ^) , both in experiment and computer simula- 
tion, is horizontal flow. While this novel effect is perhaps 
not entirely unexpected because the shape of the base 
breaks the horizontal symmetry, the complex manner in 
which the flow direction and magnitude depend on the 
many parameters that specify the system is surprising. 
The purpose of this Letter is to describe a series of simu- 
lations that resolve this apparent complexity in terms of 
the flows in horizontal strata within the vibrated layer; 
when this hitherto unreported aspect of the behavior is 
taken into account the apparently complex nature of the 
overall flow ceases to be a mystery. 

The problem of a vibrating granular layer on a flat base 
is an example of a granular system where it has proved 
possible both to validate the simulational approach based 
on suitable discrete-particle models and to use simulation 
to probe details of the dynamics that are not readily ac- 
cessible to experiment; the standing wave patterns that 
occur in this system have been studied in two [|, || and 
three |7], || dimensions. When the sinusoidally oscillat- 
ing base is given a sawtooth surface profile || horizon- 
tal flow appears; experiment and simulation both display 
the complex dependence of the flow direction and mag- 
nitude on the parameters defining the system. (While it 
is also possible to produce horizontal flow by combining 
horizontal and vertical base vibration ^ , this is distinct 
from the sawtooth problem where the vibration is entirely 
vertical.) The sawtooth base can be regarded as a kind 
of ratchet; evidence of the complexity of systems involv- 
ing potentials that produce ratchet-like effects appears in 
much simpler one-dimensional systems |n], 
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In the present series of two-dimensional simulations, 
the interaction specifying the grain shape is assumed JI^ ] 
to have a Lennard- Jones (LJ) form, with a cutoff at the 
range where the repulsive force falls to zero. For grains 



located at r, and rj this is = (48e/ry)[(cry/ri i j) 12 — 
0.5(<Jij/rij) 6 ]fij for ry < 2 1 / 6 cry and zero otherwise; 
here ry = rj — rj and cry = (<7j + Oj)/2. Ui is the ap- 
proximate diameter of grain i, although since the grains 
possess a certain degree of softness this is not precisely 
defined. The LJ interaction is strongly repulsive at small 
Tij , more so than linear overlap and Hertzian type repul- 
sions also used in granular simulation ifLsf . For conve- 
nience, we use reduced units in which length is expressed 
in terms of the diameter of the largest grains, energy 
in terms of e, and a grain with unit diameter (in re- 
duced units) has unit mass. Grain sizes are randomly 
distributed between 0.9 and 1. 

Normal and transverse viscous damping forces [ fl4[ |l5| ] 
produce the inelastic collisions and retard sliding during 
the collision. The normal force is = — 7„(ry • fij)fij. 
The transverse force is = — sign(w| J ) min(/z|/£ + 



itj\ils\vtj\)sij, where vfj = r i} 



<7jLUj)/((Ti + <Tj) is the relative tangential velocity of the 
disks, iij = z x Tij (z is the unit normal to the simulation 
plane), and is an angular velocity. The value of the 
static friction coefficient used here is /i = 0.5 and the nor- 
mal and transverse damping coefficients are j n = ~/ s = 5. 

While most granular simulations are based on particles 
with a certain amount of softness, allowing the use of 
differentiable potentials, there is an alternative approach 
involving step potentials; provided the former are not too 
soft and the flow involves some degree of fluidization the 
results of the two approaches are similar (l| [l7j. The 
previous simulations of the sawtooth problem |3|] used 
step potentials and the results were limited to thin layers; 
the softer potentials of the present work do not restrict 
the layer thickness. 

The sawtooth base is constructed from a set of grain- 
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FIG. 2: Contour plot showing overall horizontal flow veloc- 
ity as a function of / and s for h = 4 (in reduced units); 
solid and dashed curves denote positive and negative values 
respectively, and the contour interval is 0.013 (the same value 
is used subsequently). 
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FIG. 3: Flow velocity (see Fig. |) for h = 10. 

like disks positioned to produce the required profile; these 
disks oscillate vertically in unison to produce the effect of 
a sinusoidally vibrated base. The disks themselves inter- 
act with the grains using the same force laws as above; 
their diameter is 0.33 and the distance between disk cen- 
ters is half this value, producing sawtooth edges that are 
reasonably straight. The system is horizontally periodic 
and extends sufficiently far vertically to prevent grains 
reaching the upper boundary. 



Gravity also acts on the grains; here g = 5. The total 
force and torque on each grain are readily computed and 
the equations for translational and rotational motion are 
numerically integrated using the leapfrog method with 
a time step of 0.005. These and other methodological 
issues are described elsewhere Q. Each run extends 
over 1000 base cycles and the results are grouped into 10 
non-overlapping blocks to provide error estimates. 

The results described here are for a sawtooth height of 
2 and a vibration amplitude of 1 . The sawtooth shape is 
strongly asymmetric, with a gently sloping left edge and 
a very steep right edge; projected onto the base the edge 
lengths are in the ratio 99:1. The adjustable parameters 
are the base vibration frequency / and the number of 
sawteeth s (s is specified rather than the tooth width 
since it must be an integer). The initial state consists of 
grains positioned on a rectangular grid with unit spacing; 
the horizontal grid size is fixed at 90 (this is also the 
system width), and the adjustable vertical grid size is 
taken to be the nominal thickness h of the granular layer. 

Fig. |] shows a screen image (with limited visual resolu- 
tion) of a portion of a typical system for the case s = 40. 
Since the base is near its lowest point the surface wave- 
form is absent. The arrows indicate the stratified flow 
discussed below. 

Fig. H shows a contour plot of the average horizontal 
flow velocity as a function of / (frequency) and s (number 
of teeth) for layer thickness h = 4. The flow is negative 
at low s, and for certain / also at high s (correspond- 
ing to narrow gaps between teeth with space for just a 
single grain); horizontal cuts through the plot for cer- 
tain s values show velocity reversals as a function of /, 
likewise vertical cuts as a function of s. Fig. || is sim- 
ilar but for h = 10; the flow is stronger at low / and 
there is a larger negative region at high s. Examination 
of the entire range of h values considered (between 1 and 
16) reveals smoothly shifting domain boundaries separat- 
ing positive and negative flows; merely examining a few 
cuts through these two-dimensional plots will miss this 
gradual variation. These results include the complex pa- 
rameter dependence observed previously (|, H) , but with 
reversed signs because of the opposite sawtooth asymme- 
try. 

The results shown in Fig. [| deal with another two- 
dimensional cross-section through the multidimensional 
parameter space and reveal the flow dependence on s 
and h for / = 0.5. Flow reversals are seen to occur both 
as functions of s and h (including monolayers for which 
h = 1). Other values of / over the range examined (0.3 
to 0.7) again show a gradual shifting of the boundaries 
between regions of oppositely directed flow. 

Space precludes discussion of the corresponding results 
for other parameter choices (such as sawtooth height and 
asymmetry) which are similar in essence but differ in de- 
tail. A particularly drastic change is to set the tangential 
damping to zero, thereby removing the rotational com- 
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FIG. 4: Flow velocity as a function of s and h for / = 0.5. 

ponent of the dynamics. Once again the quantitative 
details change, but the overall behavior is not strongly 
affected; the main effect of rotation and the associated 
damping is to increase energy dissipation. 

Further progress in understanding the behavior re- 
quires looking beyond the overall flow and taking ad- 
vantage of the simulations to examine the dynamics in 
greater detail. An analogy can be made with thermal 
convection where the mean flow is zero, but this does 
not preclude the presence of highly structured convection 
rolls. Here, too, there are more localized flow structures 
that can only be detected by probing the spatial vari- 
ation of the flow; Fig. [j] already hinted at the kind of 
behavior that can be observed if flows in distinct hori- 
zontal strata are measured. The experiments carried out 
on this system ^ || examined just the overall flow by 
measuring tracer particle motion, and the accompanying 
simulations were confined to thin layers only, so there was 
little opportunity for observing any spatial variation. 

Studying stratified flow requires assigning grains to 
horizontal levels while allowing for the fact that grains 
can migrate vertically. Since this procedure is not 
uniquely defined, the approach adopted here is to re- 
assign the grains once every vibration cycle when the 
base reaches its lowest point. This has certain short- 
comings, particularly for higher /, because grains in the 
peaks of the surface waves are assigned to higher levels; 
however, apart from some smearing of the measurements 
between levels, no spurious effects are expected from this 
approach. 

Fig. U shows the height (z) dependence of the stratified 
flow velocity for h = 8, for several sets of / and s values; 
the thickness of each level is 1.2, which is slightly larger 
than the grain size. Fig. o shows the corresponding re- 



FIG. 5: Stratified flow velocity (v) as a function of height (z) 
above the base for h = 8; typical error bars are shown. 
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FIG. 6: Stratified flow velocity (see Fig. |) for h = 12. 

suits for h = 12. A typical velocity profile starts with 
negative or near-zero flow at the bottom level; the flow 
initially increases with z, reaches a positive maximum, 
and then starts to drop, in most cases becoming negative 
again. The curves for s = 10 (wide teeth) and s = 40 
(narrow teeth) are the most negative at the upper levels, 
while intermediate curves show strong positive maxima 
near z = 4. 

The preferred flow directions at different levels reflect 
the key features of the system. The behavior of the grains 
near the base is reminiscent of a thin layer that can flow in 
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either direction, depending on the prevailing conditions. 
On the other hand, the asymmetry of the sawteeth used 
here, which in themselves would be more likely to reflect 
a single falling grain in the negative direction, indirectly 
influences the behavior in the upper levels; this effect 
is transmitted through the intervening material which 
could well be moving in the opposite direction. When 
there is positive flow at the lower levels the competi- 
tion between these opposing effects produces the coun- 
terflowing velocity profiles. This idealized scheme is of 
course complicated by other aspects including grain ro- 
tation and the fact that the grains at the bottom can be 
permanently trapped between the teeth thereby reducing 
the effective sawtooth height. 

The explanation for the puzzling reversals in overall 
flow direction is embodied in curves such as those shown 
in Figs. H and |(| The overall flow is the density- weighted 
mean of the stratified flow velocity; its sign depends on 
the form of the velocity profile, itself a consequence of 
the relative strengths of the positive and negative flows 
occurring at different levels. 

The more familiar role of simulation, particularly as it 
applies to granular media, is in attempting to reproduce 
experimentally observed behavior. In this Letter the pro- 
cess has been inverted, and through simulation it has 
proved possible to identify a previously unknown mech- 
anism in the model granular system. On the assumption 
that the model correctly describes the relevant features of 
real granular materials this amounts to a prediction that 
calls for experimental verification. Should this kind of 
behavior be found experimentally, stratified flow would 
be an interesting addition to the extensive repertoire of 
exotic properties exhibited by granular matter. 
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